Scalar mesons in rf — > rj^n decay 
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The decay rf — > t]tt iy is studied in the framework of isobar model. It is shown, that good agree- 
ment with the experiment is achieved if ao- and cr-meson contributions are taken into account. The 
contribution of ao-meson is dominant, but cr-meson is necessary to reproduce the form of experimen- 
tal Dalitz plot. Instead of usual Breit-Wigner form of cr-meson propagator we use parametrization 
of 7T7r-amplitude, that satisfies analiticity, crossing, unitarity and chirality constraints. This ampli- 
tude has a pole in complex plane, that corresponds to cr-meson and describe experimental data on 
■ 7T7r-scattering in K e 4 decay. 
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I. INTRODUCTION 

Precision measurement of the decay 
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rf - 7?7rV (1) 

starts with the work where parameters of the matrix element of this decay in linear approximation were first 
determined with high accuracy. In this work squared matrix element |.M| 2 is expanded in Dalitz variables 

o : 

where T n o , T^o , T v are kinetic energies of n°- and rj- mesons in 7y'-meson rest frame (T^o > T n o ) and Q = T n o +T^o +T V = 
m v i — m ri — 2m n o. Now parametrization with quadratic terms proposed in Q is widely used: 

(N 1 \M\ 2 cx 1 + aY + bY 2 +dX 2 , (3) 

o . 

O^n where a, b, c, d are real numbers. 

This approach was used in subsequent works studying the decay ((T|) in neutral [|[ and charged [H modes. In both 
^ . cases one has good agreement in coefficients of expansion ([3]) and non-vanishing quadratic terms. 

• The form of the expansion ([3]) is motivated by chiral perturbation theory (ChPT). In the framework of this model 
the process r( — > 7]nir was studied theoretically in numerous works [1, H, 0, H, Q ■ It turns out, however, that leading 
order ChPT calculations gives small value of partial width of this decay and overestimate the slope parameters. The 
reason for this discrepancy is that the mass of 7y'-meson large, so leading ChPT approximation cannot be used. Final 
state interaction, on the other hand, leads to significant contribution of scalar meson resonances ao, a, /o, etc. In 
ChPT these contributions are taken into account introducing additional contact terms. The resulting width, however, 
is smaller than the experimental value. In works (To| it was shown that this problem can be solved in the framework 
of isobar model with contributions of ao, er and fo mesons taken into account. This conclusion was also confirmed in 
other works @,[Ill, [13 • 

In our paper we use this approach to analyze new experimental data on rj' — > ry7r°7r° reaction [l3l |. It is well known, 
that cr-meson contribution can not be described with usual Breit-Wigner parametrization. Analyticity, unitarity and 
chiral properties of 7r°7r° scattering, on the other hand, put strong constraints on the possible form of the amplitude of 
this process. In the works [HIGH the parametrization of the amplitude, that satisfies all mentioned above properties, 
was presented. Free parameters of this parametrization are fixed from the values of 5- wave nir scattering length and 
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FIG. 1: Typical diagrams for the process rf —* r/nn 



experimental data taken from K e ^ decay (i.e. K — > ivmr). This is the main difference of our approach in comparison 
with the other works, where simple Breit-Wigner form of er-meson propagator was used. Since we use complete 
amplitude of 7T7T-scattering in our kinematic region, contributions from other resonances (for example, /o-mesons) are 
also included automatically. In addition, for the propagator of ao-meson we use slightly modified expression, that 
takes into account the closeness of this meson to ijtt and KK thresholds. 

The rest of the paper is organized as follows. In the next section we give parameterizations for amplitudes of 
final meson interaction in nn and irr/ channels. In section III the results of the fits of these parameterizations on 
experimental Dalitz plot of rf — > rjTrir decay are presented. Discussion of our result is given in the conclusion. 



II. MATRIX ELEMENT 



In the framework of isobar model the matrix element of rj' — * 7771"° 7r° decay can be described by taking into account 
the contributions of nearest scalar resonances: ao meson in 7rf7-channel and cr-meson in 7T7r-channel (see fig|T] for 
typical diagrams). It should be noted, that there are also scalar /o-meson that can give contribution to 7T7r-scattering 
amplitude (diagram [T]d) . In the present article we use total amplitude for 7T7r-scattering in our kinematic region, 
extracted from experimental data. It is clear, that there are not only er-meson, but also /o-meson contributions in 
this amplitude. 

The amplitude of the process 77' — > rjirir can be written in the form 

A = A mi (si) + A V7] (s 2 ) + A™ (53) , 
where A^ n and A n7r are rescattering amplitudes in tttj- and 7T7r-channels respectively and 

81 = (P2+P3) 2 , s 2 = (pi + p 3 f , s 3 = (p 1 +p 2 ) 2 , 
S1 + S2 + S3 = M%,+M* + 2ml. 

Roy equations and chiral perturbation theory allow one to describe accurately the amplitude of 7T7t scattering in low 
energy region, that is allowed in considered here decay. This approach makes it possible to control the analytical 
continuation of the amplitude to the complex plane, where the pole interpreted as cr-meson is observed (l6j . 

Due to unitarity the amplitude of 7T7r-scattering should obey rather strong conditions. First of all, for s < (277V ) 
this amplitude should be real. For s > (2m 7r ) 2 up to i^if-production threshold the imaginary part of this amplitude 
should be 

In (_!_) „ fTM. 

\A n „(s)J V s 

Chiral perturbation theory, in addition, that at s = sa = the amplitude equals to zero (Adler "self-consistency" 

conditions |l7|). 

The amplitude that satisfied listed above conditions can be expressed through the series over the variable 

Vs- \J im 2 K - s 
W(S) = VW4<-* ' 

that transforms couples s-plane with cuts s < and s > (2Mk) 2 into a disc \w\ < 1 in complex w-plane (it is clear 
that w (4Af^-) = 1, w(0) = —1). The introduction of a new variable improves the convergence of the series in the 
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considered variable domain. The amplitude of 7T7r-scattering can be written in the form 

-l 



S - SA 



ID 



B Q + Biw (s) + 



where n is unknown constant that can be determined from the fit of experimental data. 

The analysis of NA48/2 results in K e 4 decay [l3| tells us, that in the presented above series one can leave only first 
two terms with the coefficients 

B = 7.4, B x = -15.1 
These values correspond to following position of cr-meson pole: 

y/s = (459 + 259i) MeV 

Let us now proceed to channel. In our kinematic region the main contribution in this channel comes from virtual 
ao-meson. It is well known that its mass is close to KK production threshold. As a result the propagator of this 
meson is different from usual Breit-Wigner form: one should introduce the self-energy corrections caused by 7r?7, -nrj 
and KK loops. It should be noted, that these corrections are significant at the pole. Our analysis, however, show, 
that in our kinematical region these corrections are small and are above current experiments accuracy. The form of 
aoirr] and aoirq' vertices, on the other hand, is important. 

These vertices can be written in several forms. First of all one can use simple point-like representation with effective 
constants g wr) and g nri i . The first constant can be determined from experimental value of ao-meson width: 

q 2 

T (ao A) = |P| « T ao = 50 -h lOO MeV, 



1.95 GeV < g nn < 2.75 GeV. 

This value agrees with g^^ = (2.46 ± 0.08 ± 0.11) GeV, presented in recent experimental work [l8j The constant g^ v i 
can be determined either from SU (3)-symmetry relations [19j or directly from fit of the considered in this article 
process. It should be noted, that SU(3) relations require information on quark structure of ao-meson, that is widely 
discussed question. For this reason we use experiment for determining this constant. 

Another form of interaction vertex is motivated by chiral perturbation theory. According to it the vertex should be 
equal to zero in the limit p v — > 0. In this case the vertex can be written in the form (p-nPi-i) l^-q for ao — * nr] interaction 
and (pnPr)') 7ir?j' for v' ~ * a o n interaction. This form was used in paper It seems more attractive, since in hit 
scattering amplitude chirality conditions are taken into account. The constant 7^ can be determined from ao — > 7777 
decay width 

5.7GCV" 1 < < 8.1 GeV -1 , 

while for determination of 7™/ constant one can use SU (3)-symmetry or distribution of 77' — > rjirir decay over Dalitz 
region. 

Thus, we will use following two variants of 7r?7-scattering amplitude: 

Ar„( S ) = TT-V< \ ^ 

s — m~ + U (s)m a 

or 

-4 wr; (s) = 77rj)'77rr; 5 j . w r • (6) 

s — m% + iL (s)m a 
III. FIT RESULTS 

Our parametrization has following unknown constants: coupling constant of a- meson with f/77 pair interaction k 
(see formula (|U)) and constants of ao-meson interaction with rjir and rfir states (g vv , g n r]' or 7^^, 7^^/ depending on 
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FIG. 2: (a) — experimental Dalitz distribution for the decay {TJ from Q; (b) — MC results with only er-meson taken into 
account; (c) — MC results with only ao-meson taken into account; (d) — do — er-interference 



the form of ao vertices). These constants will be determined using the value of rf — > ?77r 7r decay width and fit of the 
Dalitz plot of this decay, obtained in GAMS-47T experiment [3|. 

Let us first consider the point-like interaction of ao-meson ([5]). Best agreement with experimental data (with 
confidence level CL=0.52) is observed when coupling constants are equal to 

k = -4.0 ±0.3, 
9n V 9^' = (0.93 ± 0.3) GeV 2 . 

Above it was mentioned that coupling constant g nri is connected with ao-meson decay width (2 GeV < < 3 GeV) . 
As a result we get following restrictions on g wr) i coupling constant: 

0.36 GeV < gun' < 0.51 GeV. 

If we assume that ao-meson is build from light quarks only (ao ~ uu + dd) , these constants should be connected from 
77 — 77' mixing. For example, in quark mixing scheme [19j | with mixing angle $ ~ 40° the ratio of these constants 
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FIG. 3: Real and imaginary parts of the amplitudes of cr-meson (thin line) and ao-meson (dashed region) 

should be 

^ = tan$ = 0.8 

It can be easily seen, that this relation does not hold for presented above values of coupling constants g^ n and g nrl i . 
This fact is not surprising, on the other hand, since current information on internal structure of ao-meson is rather 
poor. It is, for example possible, that there is noticeable ss component in this meson, as it was mentioned in [l8j |. 

The situation is different when one use chiral form of ao vertices (expression ©)• In this case best agreement with 
experiment (CL=0.92) is observed at 

k = -4.0 ±0.3 
7^7*,' - (35±4)GeV- 2 . 

Form experimental width of ao-meson one can determine the value of tv^- 

constant (5.7 GeV 1 < j„ v < 8.1 GeV 1 ). 

It is easy to obtain the following restrictions on coupling constant 'y^t '. 

4.6 GeV" 1 < 77rt/ < 6.6 GeV" 1 . 

The relation caused by SU (3)-symmetry holds for these constants significantly better, then in the case of point- like 
a 7r?7 interaction. We would like to note, that these values are close enough (up to a sign) to results of the work 
where 7^ = 6.8 ~ 7.2 GeV -1 , = 7.4-^8 GeV -1 . It is clear, that difference in sign is insignificant, since one freely 
change the total sign of the amplitude. 



IV. DISCUSSION 



In fig[2]we present experimental Dalitz plot of 77' — > r)ir°ir° decay (figHK)j and the results of Monte-Carlo modeling 
with only ao-resonance contribution taken into account (figO)), only cr-meson contribution taken into account (fig|2fc), 
and only interference between ao and a (fig. - From these figures it is clearly seen, that the ao-meson gives the main 
contribution. The contribution of a-meson is significantly smaller, while the role of ao — cr interference is comparable 
with that of ao-meson. 

It is interesting to note, that, in spite of dominant role of ao resonance, V-distribution generated by ao-resonance 
only have opposite slope in comparison with experimental data. After inclusion of cr-meson the agreement with the 
experiment is restored. So we can conclude, that, though the contribution of cr-meson to partial width of the decay 
rj' — > 7/7T7r is small, it plays a crucial role in description of experimental data. 

This point can be clearly seen from studying contributions of different resonances to real and imaginary parts of 
the amplitude. In fig [3] we show real and imaginary parts of the amplitudes of ao and a exchange versus Dalitz 
variable Y. In our parametrization cr-meson amplitude is X independent, what can be seen from thin lines on these 
distributions. Amplitude of ao-exchange, on the other hand, depends on X, so one can see a number of lines (each 
individual line corresponds to different X values). Amplitudes shown in fig [3] agrees well with fig.ll, presented in 
paper 2]. Difference is caused by difference in parameterizations of resonance amplitudes, although we use slightly 
different parameterizations for 7T7T- and Trry-rescattering amplitudes. From this figure it is clear, that the contribution 
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FIG. 4: Squared matrix element as a function of Dalitz variable Y. Plots (from left to right) correspond to ao-meson only, 
cr-meson only, and total matrix elemet 

of ao-meson exchange dominates in the real part of the amplitude, while in the imaginary part main contribution 
comes from cr-meson exchange. 

In order to demonstrate crucial role of a- meson for description of Dalitz distribution it is useful to trace contributions 
of different mesons with presented above parameters into squared matrix element. In figure [4^, we show ^-distribution 
with only ao-meson exchange taken into account (each line in this graph corresponds to specific X value), in figQJ) 
^-distribution with only cr-meson taken into account (this amplitude does not depend on X, so only one thin line is 
present), while in fig(3t all terms of the amplitude are used. It can be clearly seen from this figure, that experimental 
slope of the ^-distribution can be achieved only if cr-meson is included. We think that this proves firmly the necessity 
of u-meson inclusion for analysis of the considered process rj — > rjTr°ir . 

V. CONCLUSION 

In our work we analyze experimental data of the reaction rf — > rjir°7r° in the framework of isobar model taking 
into account contributions of scalar mesons. We show, that it is sufficient to include ao- and er-mesons to describe 
experimental data accurately. 

In the ao-meson case one can use simple Breit-Wigner parametrization with energy-dependent width. We consider 
two different forms of ao — > tjtt and ao — > rj'n vertices and compare corresponding coupling constants with relations 
caused by SU (3)-symmetry. Our analysis show that in the case of chiral-type interaction coupling constants agrees 
S'?7(3)-symmetry much better that point-like couplings. 

The situation is more complicated for cr-meson, since the width of this meson is comparable with its mass, and 
simple Breit-Wigner parametrization cannot be used. For this reason a more accurate description of cr-meson exchange 
amplitude is required. This amplitude should satisfy a number of conditions: analiticity, unitarity, and crossing 
symmetry. In our article we use 7T7r-amplitude, obtained from fit of the experimental data of A" e 4-decay [l3j |. It 
contains a pole in complex plane at y/s = (459 + 259z) MeV, that is associated with cr-meson. 

We show, that, in agreement with previous work [2], 7r°?7-rescattering via virtual ao-meson gives main contribution 
to partial width of the 77' — ► r]Tr°ir decay. For description of the Dalitz distribution, on the other hand, it is necessary 
to take into account cr-meson contribution and its interference with ao one. If these effects are neglected, the slope 
of the Dalitz distribution in the variable Y is opposite to experimental value. Inclusion of the cr-meson corrects the 
situation. It should be noted, that cr-meson is very exotic particle, that has width comparable with its mass. For this 
reason direct experimental observation of this particle is very problematic. We think that presented in this article 
analysis of rf — ^ ry7r°7r° decay gives additional argument in favor of cr-meson existence. 
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